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We consider in this work continuous gravitational wave (GW) emission Irom non-axisymmetric radio 
pulsars. We treat in some detail the observational issues related to the known radio pulsar sample 
with the aim of unveiling the actual number of sources contributing to GW, which are likely to 
be the main contributors of GWs. It is shown that the operation of spheroidal GW detectors and 
full-size interferometers could detect this component of the radiation or impose useful limits on the 
effective oblateness of young radio pulsars. 
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I. INTRODUCTION 

Considerable experimental and theoretical efforts are 
being devoted by several groups to the physics of gravita- 
tional waves (hereafter GW). After the celebrated birth 
of neutrino astronomy [1] , GW detection may open a new 
observational window related to a long-sought prediction 
of general relativity. 

Presently, not only several resonant bars are being op- 
erated and continuously improved [2], but also the first 
full-scale wideband interferometers have been started and 
should be operational before the turn of the century [3] . 
The expected sensitivity of the advanced versions of these 
detectors in terms of the gravitational strain is of the or- 
der of /i « 10~^^ for short-lived, impulsive bursts, and 
h ~ 10~^^ for periodic ("long-lived") sources, for which 
integration times of about 10^ s are possible. In addition, 
spheroidal antennas of restricted bandwidth have been 
put forward again to complement the former, specially 
in the kHz region where they are expected to have their 
best performances [4]. This is precisely a range where 
we expect to find most of the nearest sources, namely 
rotating neutron stars in the Galaxy. The first positive 
detection(s) of GW by any of these devices will be a pro- 
found and significant achievement in physics. 

Clearly, a major concern for the above detectors is to 
understand and minimize the sources of noise which de- 
fine the sensitivity of the antennas at a given frequency. 
Furthermore, the important question of identifying a 
standard source has been put forward [5]. Among the 
possible candidates, it is expected that pulsars, which are 
known to be abundant in the Galaxy, could rank among 
the most conspicuous emitters of GW. Pulsars could have 
a time-varying quadrupole moment (and hence radiate 
GW) by either having a slight asymmetry in the equato- 
rial plane (assumed to be orthogonal to the rotation axis) 
or a misalignment between the symmetry and total angu- 
lar momentum axes, which produces a wobble in the star 
motion. In the former case the GW frequency is equal to 



the rotation frequency, whereas in the latter two modes 
are possible: one in which the GWs have the same fre- 
quency as rotation, and another in which the GWs have 
twice the rotation frequency (the first mode dominates by 
far at small wobble angles while the importance of the 
second increases for larger values). These mechanisms 
have been considered in the literature [6-14] and there 
seems to be a certain consensus for the likelihood of a 
positive detection once the antennas become operative. 

Generally speaking, the calculations and discussions 
have been focused on the detection of individual sources 
like the Tokyo experiment, devised to detect GW from 
the Crab pulsar [15]. However, very recent work along 
these lines pointed out that the radiation from the pul- 
sar population as a whole could be detectable and paid 
particular attention to the expected time modulation in 
full-scale interferometers [16,17]. On the other hand, om- 
nidirectional detectors (which may well be operative on 
shorter time scales) would not experience (by definition) 
any modulation and will be designed to have a very good 
sensitivity for continuous radiation and integration times 
of the order of several months. Thus, we have in princi- 
ple a physical component which is different from simple 
noise background. We address in this work the expected 
features of this signal, with particular emphasis on the 
difficulties created by the use of the observed distribution 
of the sources, which is affected by several observational 
biases. The goals of this first approach will be to unveil 
as much as possible the features of the true population 
contributing to this emission and to discuss what can be 
learned from the observations of the latter. 



II. PULSAR STATISTICS 

In its simplest form, the evolution of a single pulsar 
angular frequency may be written as 
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where lUq is the frequency at birth, is the character- 
istic e-folding time for the pulsar deceleration and m is 
the index of the torque power-law. From this equation, 
the so-called braking index is defined as n = ljuj/lu^ = 
1 + 2/m. As is well known, the case of pure dipo- 
lar magnetic braking corresponds to to = 1 and hence 
n = 3. On the other hand, the four presently determined 
braking indexes differ from the canonical value (2.51 for 
PSR0531 + 21, 2.24 for PSR05A0 - 69 and 2.837 for 
PSR15Q9 — 58 where the errors affect the last significa- 
tive digit, see [18] and references therein); most notably 
the recently announced value 1.4=f0.2 for the Vela pulsar 
by Lyne et al. [19]. The true reasons for these discrep- 
ancies are not clear as yet, but they may reflect substan- 
tial non-dipolar components of the magnetic field B [20], 
time- varying physical parameters like the amplitude of B 
[21], the moment of inertia / [19] or the angle between 
the magnetic and rotation axes [22,23]. For the moment, 
we shall leave m unspecified so as to allow realistic non- 
canonical braking scenarios as well. The braking will be 
determined as a byproduct of the analysis of the pulsar 
period distribution, and therefore should be considered 
in a statistical sense only and not aiming to represent 
any particular object. 

To address the expected GW signal from the pulsar 
population in the galactic disk, we should first express 
the number of observed pulsars per unit volume and pe- 
rio d interval (fN/dPdV function of puLsar-rolatcd 
parameters and observed quantities. This number can be 
written as 
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where is the; 



pulsar birthrate and Vd is the disk 
volume. The factor fg accounts for the beaming effect 
and for most of the emission models is in the range 0.1 
to 0.25. The quantity p{R, Z) is the probability to find 
a pulsar at a galactocentric distance R and at a height 
Z from the galactic plane. The pulsar observations at 
radio frequencies, where they have been widely studied 
(around 400 MHz), may be used as a "window" defining 
the "observable" volume, since we expect that detection 
decreases considerably for flux densities ^400 below 10 
mJy. As usual, the flux density (at 400 MHz) is defined 
by 5* = Lp/r"^, with Lp being the pulsar luminosity and 
r the pulsar distance to the observer. Since the radio 
sensitivity is necessarily limited by instrumental and sig- 
nal analysis [24] to those objects emitting stronger than 
Smin — lOmJy (meaning that we have an incomplete 
knowledge of the true sample, see for example [25]) the 
integration over the volume is converted into 
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To proceed we need to evaluate the integrals in cq.(3). 
To this purpose we have checked that the observed pulsar 
luminosity function Lp can be represented to a very good 
accuracy by 



AP" exp{-(3P) mJykpc^ 



(4) 

0.982 



where A = 6A x lO"^ , a = 0.37 and (3 
are the adjusted parameters. 

Let us, for the moment, assume an uniform space dis- 
tribution through the disk. Therefore using eqs.(2),(3) 
and (4) we obtain after integrating 
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To extract the value of the unknown parameters, we 
proceed to plot log dN/ dP exp{pP) as a function of log P 
for the actual observed distribution (see, for example, 
[26,27]). This yields 



log 



dN 
dP' 



expiPP) = 3.00 + OAllogP 



(6) 



with a correlation coefficient of 0.904 . Therefore, by 
comparing eqs.(5) and (6) we get m — 1.212 , and thus 
a (statistical) braking index n = 2.6, not far from the 
mean of value the directly observed braking indexes [18]. 

In order to estimate the pulsar birth rate [26], we first 
note that the constant in eq. (6) is defined by the relation 
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Prom eq.(6) one obtains 
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Performing the average of the above equation, using 
the data from reference [27] yields 
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and inserting this value in eq. (7) , one gets fs Vp = 
10^'^ yr^^, which corresponds to a pulsar birth rate of 
once every 100 to 200 yr according to the adopted value 
for the beaming factor. These results are in good agree- 
ment with the calculations performed by Lorimer et al. 
[28]. On the other hand, assuming an average initial pe- 
riod of pulsars of 10 ms, we conclude from eq.(9) that 
the average braking time is « 1022 yr and that their 
average lifetime is about 42.3 Myr. 

Now we have expressed the observed number of pulsars 
as in eq.(5), we can state that the actual number is simply 
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dN/dP times fg^, or going back to eq.(2) and changing 
, for convenience to the variable u> 
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where we have introduced A 
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With these results we turn now to the issue of calculating 
the GW emission expected in the detectors. 



factors of 0(1). In this case we obtain, after integrating 
eq.(ll) over the range 10-1000 Hz, 



< /z2>V2= 3.3x10- WJis- 



(13) 



where 0,max is the upper frequency bound of the an- 
tenna. 



IV. DISCUSSION 



III. EXPECTED AMPLITUDE OF THE GW 
EMISSION 

A population of pulsars having a disk distribution, and 
for which a slight equatorial asymmetry in each object is 
expected to be present, produce a non- vanishing value of 
the gravitational strain h?, whose average value is given 

by 

<u'->^ j,iv j,ml{£L) (II) 

where hi = ^ (§)^ £^ ^ is the angle-averaged 
contributed amplitude from each source [7,9] adequate 
for omnidirectional antennas, (sec [16,17] for a discus- 
sion of the interferometer case where the orientation is 
important), Izz is the principal moment of inertia about 
the z axis (assumed to be the same for all pulsars) and 
e" = (a — b)/{ab)^/^ is a mean over the individual ellip- 
ticities in the equatorial plane of the stars (a, b being the 
radii along the x and y axes respectively) . il = 2ttv is the 
GW frequency and is equal to twice the pulsar rotation 
frequency w and r is the distance to the pulsar. Note that 
this emission is assumed to arise from non-axisymmetric 
bodies and not from precessing ones, which have been 
addressed in dc Araiijo et al.(1994) [14] and which arc 
likely to appear as impulsive sources (in spite of several 
remarks along the years some confusion remains in the 
literature about these cases). 

Inserting eq.(lO) into eq.(ll) and integrating over the 
volume and bandwidth AO of a resonant detector we get 

< e >l/2 = 3.27 X 10-24 O0-675 AnV2 g (12) 

where fie is the central frequency of the antenna. In 
order to obtain the numerical factor, we have adopted 
a disk radius of 12 kpc and that the nearest pulsar is 
about 500 pc. Increasing or decreasing the latter value 
by a fac:tor of two produces a decrease or an increase by 
20% respectively in the coefficient of eq.(12). 

The case of broad-band interferometers is a bit more 
complicated since the signal does depend on the rela- 
tive orientation an thus the signal will be modulated (see 
[16,17]. Although we have not attempted a detailed cal- 
culation of this modulation, our model is nevertheless 
useful to give an estimate of the emission strength modulo 



From the results of the former sections we can out- 
line an optimal strategy to search for this signal. Let 
us consider first the case of a pair of resonant spheroidal 
antennae [29]. We assume a design capable of cover- 
ing the frequency range 800-1000 Hz and a strain noise 
^JSh 10~^^ Hz~^^'^, which seems reasonable for the near 
future. In fact, it would be possible to go as low as 600 
Hz for the central frequency if the spheroidal detector 
can be made of niobium instead of a cheaper material 
like Al alloys. Such a reduction is important when the 
whole set of possible sources is considered and is under 
study by some experimental groups [29]. 

For each of these omnidirectional detectors, the posi- 
tive identification of this GW radiation at, say, 4cr level 
requires a minimum detectable amplitude 

C 1/2 

hc^H^) (14) 

where we shall assume that the strain noise ^/Sh is 
dominated by thermomechanical forces [4,13] and r is 
the integration time ~ 10^ s. Imposing reasonable tech- 
nological improvements for a pair of Al alloy " 4*'* " gen- 
eration antenna, the minimum detectable amplitude is 
he > IQ-^^. Comparing eq.(12) and eq.(14) we deduce 
that the emission from the ensamble would be detectable 
by the array if the (mean) characteristic ellipticity of pul- 
sars satisfies 

£ > 3.6 X 10-^ (15) 

That is, the array would be sensitive to the pulsar pop- 
ulation rotating faster than P ~ 2.4 ms. These figure 
change toe > 4.7x10-'' and P simeq3.6m,s for a nio- 
bium array with a central frequency of 600 Hz, which 
gives an idea of the range to be explored by these de- 
vices. 

Since we have always referred to the radio pulsar pop- 
ulation, for consistency of the picture we shall demand 
that < Izz^^ , where the right term represents the 
energy rate being extracted from the rotating neutron 
star. Using the well-known expression for Eg^ derived 
by Ferrari and Ruffini [6] 

Eg^ = ^^^lle'u:' (16) 

we get that the GW emission is not dominant provided 
that [13] 
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£ < 1.9 X 10^ (PP3)^/^ (17) 

For those target pulsars contributing to the pair of 
spheroidal antennae we get from eqs.(l) and (8), an av- 
erage initial deceleration rate is P w 5 x 10""'^'^, and 
thus their average equatorial deformation must satisfy 
e < 1.5 X 10~^ in order to brake by electromagnetic 
emission. This upper limit is a factor of 40 higher then 
the lower limit implied by eq.(15) an)d, in a sense, the 
window to be probed by the detectors. It is interesting 
to note that the window does not shrink (in fact becomes 
slightly larger) if one considers the lower frequency nio- 
bium array. Although we shall not be concerned with 
the precise origin of the ellipticity, we note that if the 
crust of neutron stars is solid, then its shape may not 
necessarily be axisymmetric under the effect of rotation, 
as it would be the case for a fluid. More recent stud- 
ies suggest strongly that most of the neutron star is in 
a liquid phase [30], and in this case other mechanisms 
must be invoked to produce triaxiality. Several other 
mechanisms, however, may induce non-axisymmetric de- 
formations above the detectability bound of eq.(15) as 
discussed in refs. [31-34] (in particular, the presence of 
stochastic magnetic fields in the mantle of the neutron 
star or the existence of a type II superconductor [35]). 

The possible existence of a subpopulation of silent ra- 
dio pulsars whose spin-down is driven by GW [13] is also 
interesting and would add to ft., although we have not 
based any of our predictions on them. We also remark 
that, given the uncertainties governing our knowledge of 
the young pulsar sample, it is not guaranteed that the ini- 
tial periods of the newborn pulsars are as short as 2.4 ms 
(or 3.6 ms in the case of a niobium array). In fact, a direct 
backward extrapolation of the present Crab period gives 
an initial value of 4 — 5 times of that, or about 10 ms. If 
this is the generic case it is unlikely that any omnidirec- 
tional array could be constructed to detect these sources, 
although interferometers will certainly do due to their 
braod-band coverage and maximum sensitivity at lower 
frequencies. A complicated temporal modulation pattern 
will arise in the latter, which may nevertheless help to ex- 
tract this signature (sec [16]). An estimate of the signal 
strength can be obtained from our model by comparing 
eq.(13) with the maximum sensitivity region of the ad- 
vanced interferometers, giving the detectability condition 

£ > 7 X 10"^ (18) 

to be compared with the upper bound s < 8.5 x 10~^. 
The region to be explored is now about two orders-of- 
magnitude in the £ parameter. 

To understand what do the detectability conditions re- 
ally mean we have collected in Table 1 the relevant num- 
bers of five "top candidates" according to the standard 
procedure of assuming that all the observed P is due 
to GW emission. This bold assumption yields a max- 
imum ellipticity Em and therefore the maximum ampli- 
tude hm to be expected for that individual sources. Table 



1 stresses the importance of exploring the hidden sources 
we have been referring to: while the largest hm belong to 
young pulsars like Crab and Vela, the required £„ seem 
to be rather extreme (a hint for the incorrectness of the 
"all GW" assumption); on the other hand ms pulsars like 
PS'P1957 -I- 20 can not be strong GW emitters because 
their total energy loses are known to be small. There- 
fore the ideal emitters would be rapidly rotating, high 
P young pulsars which may populate the upper corner 
of the P — P diagram, having the additional advantage 
of possibly falling in both the spheroidal and interfero- 
mctric detector ranges. It can be said that the search 
for an optimal experimental situation has forced us to 
address the statistical expectations for these yet undis- 
covered sources. 



Source em hm vaw [Hz] 

Vela 1.8 X 10"^ 1.9 x 10"^^ 22.5 

Crab 7.5 x 10"'' 1.5 x 10"^** 60.6 

Gomiriga 2.3 x 10"^ 1.2 x 10"^* 8.4 

PSR1509-68 1.4 x 10"^ 5.8 x 10"^^ 13.2 

PSR1957-I-20 1.6 x 10"^ 1.7 x IQ-^^ 1244 



TABLE I. Upper limits to the GW emission from individ- 
ual sources (see text) 



4 



It is to be noted that, according to their small asym- 
metry as deduced from their observed features [13], the 
millisecond pulsar subpopulation would not contribute to 
the ensamble radiation. Consistently, this subpopulation 
has not been included when considering the statistical 
features exposed in Sec. 2. This does not mean, however, 
that ms pulsars arc not interesting for GW searches (see, 
for example, [13,14]), but rather that our calculations do 
not depend on them. We believe that the accumulated 
evidence for a distinct evolutionary path for these objects 
is enough for a separate treatment of their emission. 

As a final remark it is important to stress that, in spite 
of the statistical treatment given in this work, < h? >i/2 
will be dominated by the youngest pulsars in the Galaxy 
because of the rather strong dependence of the gravita- 
tional emission on uj. For example, from our statistical 
analysis we found that the total number of pulsars with 
P < 4s is around 10^, but only ^ 25 objects would be 
responsible for 70% of the total gravitational radiation in 
the interferometers. We would also predict in our model 
that only ~ 3 of them would be detected at radio fre- 
quencies (400 Mhz) in the period range 10-30 ms. The 
true situation may be more optimistic than we think. 
According to Camilo [36] , a clue about the young pulsar 
population may be hidden in the P — P diagram. His ar- 
gument is precisely that evolution along a B = constant 
line may require the presence of many young pulsars with 
P > 10^^'*s/s and P < 50 ms which remain undiscov- 
ered because of selection effects and are missed in the ex- 
isting searches. This independent evidence reinforces the 
motivation of our statistical approaches and calls for fu- 
ture refinements. Detection of this GW radiation would 
be an important check for the first GW observatories, and 
even a non-detection would be useful to set constraints on 
the shape an on the statistics of young galactic pulsars. 
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